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(nondiver gent)  background.  This  representation  corresponds  approximately 
to  the  use  of  the  mean  range  of  local  density  correlations  and  of  the 
macroscopic  order  as  fundamental  (parametric)  variables  to  locate  positions 
in  thermodynamic  space.  The  representation  thus  obtained  is  used  to  pre- 
dict effects  of  gravity  and  the  shape  of  contours  of  equal  wave  scattering 
close  to  the  critical  point.  The  predictions  are  investigated  experimentall; 
by  studying  the  dependence  of  the  scattering  of  laser  light  on  chemical 
potential  (height)  and  temperature.  Satisfactory  agreement  is  found  to 
within  the  precision  of  the  measurements.  In  addition,  it  is  verified 
experimentally  that  contours  of  equal  opalescence  have  the  same  shape  as 
contours  of  equal  mean  relaxation  time  of  local  density  fluctuations.  This 
demonstrates  that  the  relaxation  time  and  the  correlation  length  obey 
equivalent  scaling-law  homogeneity  relationships  close  to  the  critical 
point.  Some  discrepancies  among  existing  experiments  and  between  experi- 
ment and  theory  are  noted  and  discussed  briefly.  In  particular,  there 
remain  unexplained  discrepancies  between  the  values  of  the  critical-point 
exponents  found  here  and  those  predicted  by  the  Ising  model. 
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The  object  of  this  work  was  to  use  the  methods  of  elastic 
and  inelastic  light  scattering  to  study  critical  point  fluctua- 
tions and  to  carry  through  theoretical  analysis  to  embrace  the 
new  experimental  results  in  an  improved  theory  of  critical  point 
phenomena.  The  principal  results  of  this  work  have  been  published 
in  The  Fhyslcal  Review  (Phys.  Rev.  A 11,  1706-1723  (1975),  article 


by  J.A.  WyLte  and  B.S.  Maccabee) . A reprint  of  that  article  is 
attached  to  this  report.  Other  publications  that  have  resulted 
from  this  work  or  are  in  preparation  are  listed  at  the  end  of  this 

repdrt. 
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A principal  focus  of  the  present  work  was  to  study  the  fluctua- 
tions in  a field.  To  this  end,  a contour  of  constant  scattered 
intensity  was  measured  near  the  critical  point  and  the  results  were 
related  to  scaling  law,  parametric  model  descriptions  of  the  critical 
region. A The  requirement  of  consistency  with  bulk  thermodynamic 
measurements  performed  farther  from  the  critical  point  led  to  a 
substantially  unique  parametrization  of  the  critical  region.  This 
result  should  be  useful  for  estimating  critical  point  behavior  to  be 
expected  still  closer  to  the  critical  point  and  for  comparison  with 
theories  of  the  critical  region  currently  being  developed.  On  the 
basis  of  the  new  measurements  and  analysis  an  experiment  is  being 
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designed  to  study  the  coherence  length  (or  second  moment)  of  the 
critical  point  fluctuations  below  Tc  and  in  a field,  partly  to 
test  ideas  published  recently  by  Tarko  and  Fisher.  (See  H.B. 

Tarko  and  M.E.  Fisher,  Phys.  Rev.  B 11,  1217  (1975)). 

A related  objective  of  the  present  work  was  to  develop  an 
empirical  description  of  the  breakdown  of  scaling  law  behavior 
farther  from  the  critical  point,  using  the  parametric  model 
established  by  the  light  scattering  experiments  performed  close 
to  the  critical  point  as  a first  approximation.  This  effort  has 
resulted  in  a non-symmetric  parametric  equation  of  state  for  the 
extended  critical  region  of  carbon  dioxide  which  describes  rather 
accurately  the  observed  susceptibility  out  to  temperatures  100° 
above  Tc  and  densities  in  the  range  0.5PC  to  2.0PC.  The  description 
thus  obtained  is  being  used  in  ongoing  work  to  facilitate  a more 
extensive  analysis  of  critical  point  anomalies  in  the  extended 
critical  region. 

The  following  publications  have  resulted  from  work  covered  by 
this  report: 
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I.  Article  in  print 

1.  THERMODYNAMIC  ANOMALIES  NEAR  THE  CRITICAL  POINT  OF  C02, 

J.A.  White  and  B.S.  Maccabee,  Phys.  Rev.  A 11,  1706- 
1723  (1975). 

II.  Brief  reports 

1.  SPINODAL  BEHAVIOR  OF  PARAMETRIC  MODELS,  J.A.  White  and 
B.S.  Maccabee,  Bull.  Am.  Phys.  Soc.  20^,  651-2  (1975). 

2.  THE  PRESSURE  ALONG  THE  CRITICAL  ISOCHORE  OF  CARBON 
DIOXIDE,  B.S.  Maccabee  and  J.A.  White,  Bull.  Am.  Phys. 

Soc.  20,  652  (1975). 

3.  THE  SPECIFIC  HEAT  OF  C02  NEAR  THE  CRITICAL  POINT,  D. 

Madison  and  J.A.  White,  Bull.  Am.  Phys.  Soc.  21^,  686 
(1976).  (The  statement  "g  = 1.31"  printed  in  the  abstract 
should  be  corrected  to  read  "g  = 1.28".) 

4.  LIGHT  SCATTERING  TEST  OF  THE  PARAMETRIC  EQUATION  FOR  C02 
VERY  NEAR  THE  CRITICAL  POINT,  P.  Kang  and  J.A.  White,  Bull. 
Am.  Phys.  Soc.  21,  685-686  (1976).  (The  statement  ”g  = 1.30” 
printed  in  the  abstract  should  be  deleted.) 


III.  In  press  and  in  preparation 

1.  A NON-SYMMETRIC  PARAMETRIC  EQUATION  OF  STATE  FOR  THE 
EXTENDED  CRITICAL  REGION  OF  CARBON  DIOXIDE,  B.S.  Maccabee 
and  J.A.  White,  Phys.  Lett.  A (in  press). 

2.  LIGHT  SCATTERING  NEAR  THE  CRITICAL  POINT  OF  CO2,  P.  Kang 
and  J.A.  White  (in  preparation;  based  on  the  1976  PhD 
dissertation  of  P.  Kang). 


■ 


I 


Supported  by  The 

contract  number  NQOoi^k  Naval  B®3earoh, 

view  a vb l r^rr 5 


PHYSICAL  REVIEW  A TO  ?99P^4Tl7.5ffi-fl,ffl71l  5 ’ MAY  1975 

<£  $; 

Thermodynamic  anomalies  near  the  critical  point  of  CCJy  'o  <% 

John  A.  White  and  Bruce  S.  Maccabee*  ^ ^ o 

> yf  A 

The  American  University,  Washington.  D.  C.  o . 

(Received  12  August  1974)  ^ 
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alescence  of  CO;  are  analyzed  to  determine  a parametric,  scaling-law  representation  for  the  critical  0 ^ 
region.  A representation  of  the  form  of  the  linear  model  of  Schofield  is  found  to  be  consistent  with  x**  © 
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= 0.09,  and  for  a slowly  varying  (nondivergent)  background  This  representation  corresponds  ap-  ® ^ 

proximately  to  the  use  of  the  mean  range  of  local  density  correlations  and  of  the  macroscopic  order  ^ 

as  fundamental  (parametric)  variables  to  locate  positions  in  thermodynamic  space.  The  representation  b 

thus  obtained  is  used  to  predict  e.fects  of  gravity  and  the  shape  of  contours  of  equal  wave  scattering 
close  to  the  critical  point.  The  predictions  are  investigated  experimentally  by  studying  the  dependence 
of  the  scattering  of  laser  light  on  chemical  potential  (height)  and  temperature  Satisfactory  agreement  A 

is  found  to  within  the  precision  of  the  measurements  In  addition,  it  is  verified  experimentally  thai  <* 

contours  of  equal  opalescence  have  the  same  shape  as  contours  of  equal  mean  relaxation  time  of  local 
density  fluctuations.  This  demonstrates  that  the  relaxation  time  and  the  correlation  length  obey  equiv- 
alent scaling-law  homogeneity  relationships  close  to  the  critical  point  Some  discrepancies  among  exist 
ing  experiments  and  between  experiment  and  theory  are  noted  and  discussed  briefly  In  particular, 
there  remain  unexplained  discrepancies  between  the  values  of  the  critical-point  exponents  found  here 
and  those  predicted  by  the  Ising  model 


I.  INTRODUCTION 
A.  Background 

As  the  temperature  of  a simple  fluid  is  in- 
creased, the  crisp  distinctions  of  macroscopic 
phase  that  separate  gas  from  liquid  at  low  tem- 
peratures fade  until,  at  the  critical  point  and  at 
higher  temperatures,  the  fluid  no  longer  exhibits 
a discontinuous  change  of  phase  with  change  of 
pressure.  It  can  then  no  longer  be  characterized 
unambiguously  as  either  liquid  or  gas. 

In  the  vicinity  of  the  critical  point,  new  features 
become  prominent.  On  a macroscopic  level,  these 
include  critical  opalescence  and  anomalies  (di- 
vergences) of  the  compressibility,  specific  heat, 
thermal  conductivity,  and  diffusivity.  Micro- 
scopically, fluctuations  become  an  increasingly 
dominant  feature.  Local  deviations  from  average 
density  (fluctuations)  are  always  present,  fleet  - 
ingly,  within  each  phase,  as  a result  of  thermal 
agitation.  However,  as  a critical  point  is  ap- 
proached the  correlations  among  local  deviations 
at  different  points  in  a fluid  increase  appreciably 
both  in  spatial  extent  and  temporal  duration.  In 
effect,  each  phase  develops  bubbles  of  above  and 
below  average  density.  These  "bubbles,”  of  ir- 
regular shape  and  diffuse  boundaries,  increase 
substantially  in  size  and  decrease  markedly  in 
mobility  as  the  critical  point  is  approached.  The 
growth  of  these  locally  correlated  deviations  is  re- 
sponsible for  anomalous  static  and  transport  prop- 
erties in  the  vicinity  of  critical  points. 


To  best  represent  the  thermodynamic  properties 
of  a fluid  in  the  vicinity  of  its  critical  point,  an 
equation  of  state  should  be  Lised  upon,  or  at  least 
compatible  with,  the  anomalous  properties  which 
arise  from  the  growth  in  spatial  extent  of  the  local- 
ly correlated  deviations.  Such  an  equation  of  state 
could  be  one  in  which  the  macroscopic  thermo- 
dynamic variables  (density,  temperature,  chemi- 
cal potential,  etc.)  are  expressed  in  terms  of  one 
parameter  which  characterizes  the  locally  corre- 
lated deviations  from  average  order  and  another 
which  characterizes  the  average  order  for  a given 
value  of  the  first  parameter.  It  is  known  that  im  - 
portant  quantitative  features  of  the  thermodynamic 
anomalies  (the  strengths  of  the  divergences)  are 
related  to  a single  quantity  that  characterizes  the 
localized  deviations  from  average  density  (a  spe- 
cific type  of  order).  This  quantity  is  the  mean 
range  of  cooperation  (correlation  length)  of  these 
deviations,  or  the  mean  “bubble”  radius.  The 
magnitude  of  the  compressibility  anomaly  and  the 
strength  and  anisotropy  of  scattering  responsible 
for  critical  opalescence1  have  been  related  to  this 
range, 2,1  as  has  also  the  strength  of  the  divergence 
of  the  thermal  diffusivity.'  If  a thermodynamic 
variable  that  is  closely  related  to  the  correlation 
length  of  density  deviations  is  chosen  as  one  of 
the  above  parameters,  then  the  other  parameter 
should  be  closely  related  to  the  average  density. 

To  incorporate  these  parameters  into  an  equa- 
tion of  state  for  the  fluid,  it  is  necessary  to  make 
use  of  a transformation  from  conventionally  mea- 
sured parameters,  such  as  temperature,  pres- 
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sure,  and  density,  to  new  variables  which  are 
related  to  the  correlation  range  and  the  average 
density.  It  will  be  convenient  for  this  purpose  to 
employ  a parametric  representation.5’8  A pair 
of  parametric  variables,  r and  &,  will  be  defined 
as  functions  of  temperature  T and  density  p.  They 
will  be  chosen  so  that  r ( the  “remoteness”  from 
the  critical  point)  is  closely  related  to  the  correla- 
tion range  of  the  fluctuations  (increasing  as  that 
range  decreases)  and  9 (“order  parameter”)  is  an 
increasing  function  of  the  average  density  at  con- 
stant r that  defines  the  macroscopic  phase  of  the 
fluid.  Together  9 and  r suffice  in  place  of  p and  T 
to  specify  any  equilibrium  thermodynamic  state. 

Several  parametric  models  have  been  proposed,1"9 
the  simplest  of  which  is  known  as  the  (unrestricted) 
linear  model.8'7  In  this  work  we  have  compared  the 
predictions  of  this  model  with  C02  thermodynamic 
and  light -scattering  data.  However,  in  doing  so, 
we  have  imposed  a new  requirement  on  the  model. 
We  have  required  that  the  parameter  r be  as  nearly 
as  possible  a simple  (exponential)  function  of  the 
susceptibility,  x = (p2Kr),  because  this  suscepti- 
bility is  closely  related  to  the  density-fluctuation 
correlation  range2,3  in  the  immediate  vicinity  of  the 
critical  point  where  the  linear  model  is  expected 
to  be  applicable.  The  consequences  of  this  require- 
ment are  as  follows:  (i)  The  constant  known  as  b2 
in  the  model  is  smaller  than  the  value  proposed  by 
Schofield,  Litster,  and  Ho  (SLH)7;  (ii)  as  a result 
of  (i)  the  model  has  a "true”  spinodal  behavior  in 
that  the  compressibility  diverges  along  a curve 
that  lies  inside  the  coexistence  curve;  (iii)  sus- 
ceptibility and  chemical -potential  data  in  the  vicin- 
ity of  the  critical  point  are  well  represented  by  the 
model;  and  (iv)  the  available  specific -heat  data 
can  also  be  well  represented  by  the  model.  Thus, 
by  using  PVT  data,  specific -heat  data,  light-scat- 
tering data  that  are  found  in  the  literature,  and 
light-scattering  data  that  are  presented  here  for 
the  first  time,  we  have  been  able  to  adjust  the  con- 
stants in  the  linear-model  equation  so  that  it  rep- 
resents the  available  data  to  within  about  a 5% 
accuracy.  We  believe  that  this  is  the  first  time 
such  good  over -all  agreement  with  critical-point 
data  has  been  achieved. 

A detailed  plan  of  this  paper  and  summary  of 
our  principal  new  results  is  presented  in  the  fol- 
lowing synopsis. 

B Synopsis 

In  the  first  of  the  sections  to  follow  (Sec.  II) 
bulk  measurements  of  pressure,  volume,  and 
temperature'0  are  analyzed  to  determine  contours 
of  constant  susceptibility."  These  contours,  to- 
gether with  the  coexistence  curve'2  and  the  re- 


quirement we  impose — that  contours  of  constant 
susceptibility  approximate  contours  of  constant 
r — suffice  to  determine  three  of  Schofield's  linear- 
model  parameters,8  62,  g,  and  ji,  to  two  significant 
figures.  Agreement  for  this  choice  of  parameters 
with  some  of  the  experimental  data  is  exhibited  in 
Table  I.  Reasons  are  stated  for  the  small  dis- 
crepancies apparent  in  this  table.  In  Appendix  C 
it  is  shown  in  detail  how  the  data  discussed  thus 
far  (susceptibility  contours  and  coexistence  curve) 
lead  naturally  to  the  particular  linear-model  pa- 
rametric equations  of  state  (summarized  in  Ap- 
pendix B)  that  are  used  in  the  rest  of  this  investi- 
gation. 

Next,  in  Sec.  Ill,  the  temperature  dependence 
of  the  compressibility  at  critical  density,  as  de- 
rived from  PVT10  and  light -scattering  data, 13,14  is 
used  to  determine  independently  of  62  and  ii  accu- 
rate numerical  values  of  the  remaining  two  model 
parameters,  a and  >.  Reasons  are  presented  to 
explain  differences  between  previously  proposed 
values  for  y.  The  value  of  a,  which  is  the  least 
precisely  determined  of  the  model  parameters,  is 
found  to  be  constrained  by  the  compressibility  data 
to  be  a certain  (empirical)  function  of  two  of  the 
other  parameters,  g and  y.  The  linear -model 
parameters  determined  thus  far  give  a predicted 
ratio  of  susceptibilities  above  and  below  Tc  that 
agrees  with  experiment.13 

Numerical  values  for  the  linear -model  param- 
eters found  in  other  recent  investigations15'18  are 
noted.  They  differ  considerably  among  themselves 
and  some  of  them  differ  considerably  from  those 
obtained  here.  Some  reasons  for  the  discrepancies 
are  stated. 

The  work  to  this  point,  described  in  Secs.  II  and 
III,  leads  to  a complete  specification  of  the  thermo- 
dynamics of  COj  in  a linear -model  approximation. 
The  adequacy  of  this  model  is  tested  in  the  follow- 
ing sections  by  comparing  in  detail  predictions  of 
the  model  (Appendix  B)  with  several  different  types 
of  experimental  data. 

First,  in  Sec.  IV,  a more  detailed  comparison 
is  made  with  the  PVT  data.10  The  model  permits 
susceptibilities10  and  chemical  potentials11  derived 
from  pressure  measurements  to  be  normalised  to 
two  universal  curves  (Figs.  4 and  5).  Agreement 
with  model  predictions  is  to  within  the  experi- 
mental uncertainty  of  ~5%  for  r<0.03  and  |0|<O.8. 
For  larger  values  of  r and  9,  the  breakdown  of 
perfect  scaling-law  symmetries  becomes  appre- 
ciable. Reasons  for  the  departures  from  perfect 
symmetry  are  discussed.  It  is  noted  that  no  error 
outside  experimental  uncertainties  is  made  for 
any  6 for  our  choice  of  linear-model  parameters, 
including  especially  our  choice  for  62,  for  small  r. 
(A  model  equation  of  state  that  incorporates  cor- 
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rection  terms18  to  describe  the  asymmetries  at 
larger  r has  been  fitted  to  the  C02  susceptibility 
data  and  been  found  to  work  out  to  r = 0.3  and  8 
= ±1.0.19) 

Model  predictions  are  compared  with  specific  - 
heat  data20  in  Sec.  V.  The  comparison  leads  to  a 
further  refinement  of  the  model  parameters.  To 
make  this  refinement,  it  is  necessary  to  assume 
that  the  exponent  a that  describes  the  divergence 
of  the  specific  heat  is  smaller  than  the  value  0.125 
originally  found20  to  fit  the  specific -heat  data.  It 
is  noted  (Appendix  D)  that  values  as  small  as  0.08 
can  result  in  a fit  to  the  measured  specific  heat 
to  within  ~1%.  Using  a=0. 09,  we  find  an  excellent 
fit  to  both  branches  of  the  specific -heat  curve 
(above  and  below  Tc)  for  choices  of  the  other  pa- 
rameters of  the  model  compatible  with  the  PVT 
and  light -scattering  data,  and  compatible  with  our 
restriction  on  the  choice  of  b2 . For  our  best  fits, 
the  background  contributions  to  the  specific  heat 
are  the  same  above  and  below  Tc.  Our  best  fits, 
however,  leave  a discrepancy  ~2%  very  close  to 
Tc  that  cannot  be  removed  by  further  small  ad- 
justments in  the  values  of  the  parameters. 

A final  choice  of  model  parameters  is  made  and 
exhibited  in  Sec.  V (Table  III).  The  values  for  the 
parameters  stated  there  are  compared  with  values 
proposed  for  the  same  parameters  in  other  recent 
investigations.15'19  The  parameter  values  found  in 
the  present  investigation  give  an  order  of  magni- 
tude improvement  in  the  fit  to  the  specific -heat 
data,  and  in  particular  to  the  ratio  of  specific 
heats  measured  above  and  below  Tc.  Agreement 
over-all,  with  specific -heat,  PVT,  and  light- 
scattering data,  is  good  enough  (~3%  with  much  of 
the  most  accurate  data)  to  justify  a confidence  that 
the  model  gives  an  essentially  correct  description 
of  all  critical -point  phenomena,  at  least  when  r 
is  sufficiently  large  to  assure  that  the  compres- 
sion of  the  fluid  under  its  own  weight  is  of  second- 
ary importance.  In  the  sections  that  follow,  the 
model  thus  specified  is  used  closer  to  the  critical 
point  (for  r<10"4)  to  interpret  experiments  per- 
formed to  study  effects  of  this  compression. 

The  variation  of  density  with  vertical  position 
in  the  fluid  is  considered  in  Sec.  VI.  Height  pro- 
files predicted  by  the  model  equations  for  iso- 
thermal conditions  are  compared  with  profiles 
determined  by  Schmidt  and  Straub21'22  from  mea- 
surements of  the  refractive  index  of  C02  as  a 
function  of  height.  There  is  serious  disagreement, 
amounting  to  a factor-of-3  discrepancy  in  the 
“best -fit”  value  to  use  for  the  model  parameter  a. 

A detailed  consideration  (Appendix  F)  shows  that 
the  discrepancy  can  be  explained  by  the  presence 
of  a thermal  gradient  in  the  experiments  of  ~0.6 
mdeg/cm  (top  of  the  cell  warmer  than  the  bottom). 
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A thermal  gradient  of  this  magnitude  seems  to 
have  been  within  the  thermal  tolerances  of  the 
experiments. 

In  Secs.  VII  and  VIII  we  turn  to  the  analysis  of 
some  light -scattering  measurements  performed 
very  close  to  the  critical  point  and  reported  here 
for  the  first  time.  These  measurements  have  been 
performed  in  a region  (rs.10"1)  where  effects  of 
the  compression  of  a fluid  under  its  own  weight 
are  prominent. 

In  Sec.  VII  we  consider  the  variation  of  the  sus- 
ceptibility of  the  fluid  at  equilibrium  with  tem- 
perature and  vertical  position.  The  variation  of 
the  susceptibility  predicted  by  the  model  equations 
is  compared  with  the  results  of  light-scattering 
measurements.23  The  intensity  of  scattering  is 
expected  to  be  proportional  to  the  susceptibility1  '2 
for  infinitely  thin  samples.  However,  for  the  0.5- 
cm -thick  samples  used  in  the  experiments,  ex- 
tinction is  large,13  with  the  result  that  measured 
intensity -height  profiles  uncorrected  for  extinction 
differ  appreciably  from  predicted  susceptibility 
profiles  (Fig.  6). 

In  Sec.  VIII,  a novel  method23  for  making  light- 
scattering measurements  close  to  the  critical 
point  is  described.  In  this  method  the  temperature 
and  vertical  position  within  the  fluid  are  both 
changed  to  map  a contour  of  constant  scattered  in- 
tensity. The  shape  of  this  contour  is  largely  in- 
sensitive to  effects  of  extinction  and  tests  directly 
the  shape  of  the  underlying  thermodynamic  surface 
without  the  need  to  make  corrections  for  extinc- 
tion. The  shapes  of  predicted  and  observed  con- 
tours are  compared  (Fig.  8),  and  agreement  is 
found  to  be  good,  to  within  the  '10*?  accuracy  of 
existing  measurements .23  Thus,  no  evidence  is 
found  to  this  accuracy  for  a breakdown  of  the  mod- 
el or  of  a need  to  use  different  parameter  values 
very  close  to  the  critical  point. 

Additional  data  are  presented  in  Sec.  VIII  on 
the  spectrum  of  the  scattered  light.  The  inverse 
of  the  width  of  the  spectrum  is  proportional  to  the 
mean  time  for  growth  and  decay  of  the  spontaneous 
fluctuations  of  density  within  the  fluid,  and  hence 
is  a measure  of  a transport  property  (diffusivity) 
of  the  fluid.  Contours  of  constant  spectrum  (con- 
stant relaxation  time)  have  been  mapped  and  found 
to  be  of  the  same  shape  as  contours  of  constant 
scattered  intensity.  Thus,  to  within  the  accuracy 
of  the  measurements  reported  here,  contours  of 
constant  transport  phenomena  (diffusivity)  coincide 
with  contours  of  constant  static  properties  (sus- 
ceptibility). These  measurements,  performed 
close  to  the  critical  point  (rs  10‘4),  contribute  to 
confidence  that  a single  critical -point  scaling-law 
theory,  approximated  for  C02  by  the  model  equa- 
tions and  choice  of  b2  employed  in  this  investiga- 
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tion,  correctly  describes  the  cooperative  phenom- 
ena responsible  for  all  anomalies,  dynamic  and 
static,  present  in  the  vicinity  of  a critical  point. 

A summary  of  the  principal  general  conclusions 
of  this  investigation  is  presented  in  Sec.  DC. 

II.  COEXISTENCE  CURVE  AND  SUSCEPTIBILITY 
CONTOURS 

Extensive  and  accurate  PVT  measurements 
made  several  decades  ago10  establish  that  the 
thermodynamic  surface  of  C02  has  the  general 
shape  shown  in  Fig.  1.  Of  particular  interest  are 
the  isoclines  of  that  surface,  which  define  contours 
of  constant  susceptibility  (see  the  summary  of 
definitions  given  in  Appendix  A).  These  contours 
when  viewed  in  the  T,  p plane  [Fig.  1(b)]  have  the 
form  of  a family  of  inverted  parabolas.  Experi- 
mentally determined  contours  are  shown  in  Fig. 

2. 11  The  data  in  that  figure  cover  a rather  wide 
range  of  temperature  (10-140°C)  and  density 
(60-460  amagat)  in  the  single-phase  region  around 
the  critical  point  (rc  = 31°C,  pc  = 237  amagat).  The 
data  cover  a large  enough  range  to  show  depar- 


FIG.  1.  Thermodynamic  surface  in  the  vicinity  of  the 
critical  point.  (a)  The  real  (physical)  surface  is  bounded 
by  the  heavy  line.  This  line  consists  of  contours  of  con- 
stant AT  and  contours  of  constant  Am  = M(P,  T)  - p(pc,  T). 
The  coexistence  curve  lies  in  the  T-p  plane.  The  cross- 
hatched  regions  lie  "below'’  the  surface.  The  lines  which 
lie  in  and  illustrate  the  curvature  of  the  surface  are  iso- 
therms. The  short-long  dashed  line  is  an  "isocline'’ — 
a line  of  constant  slope,  lap/fi p)T.  (b>  The  constant 
Op/dply  line  projected  onto  the  T-p  plane,  (c)  The 
same  isocline  viewed  along  the  temperature  axis.  Id) 

The  isocline  viewed  along  the  density  axis.  The  critical 
point  is  at  the  origin  of  each  figure;  the  figures  are  not 
drawn  to  scale. 


tures  from  perfect  symmetry  at  substantial  dis- 
tances from  the  critical  point. 

Significant  features  of  the  susceptibility  data  of 
Fig.  2 close  to  the  critical  point  are  described 
most  simply  in  terms  of  parametric  variables, 
illustrated  in  Fig.  3.  An  approximate  representa- 
tion of  these  features  is  given  by  contours  of  con- 
stant r and  constant  6 in  the  linear -model  para- 
metric representation  of  Schofield,*  illustrated 
in  Fig.  3(b)  for  particular  choices  of  the  param- 
eters of  the  model.  The  equations  of  the  contours 
shown  in  Fig.  3(b)  are 

Ap  =gpcrBe,  (la) 

AT  = Tcr(l  - 6*0*).  (lb) 

Here  pc  and  Tc  are,  respectively,  the  density  and 
temperature  at  the  critical  point,  Ap  -p  -pc, 

A T -T  - Tc,  and  g,  ii,  and  are  three  dimension- 
less constants,  withb2>l.  The  dimensionless  con 
stants  determine  the  shape  of  the  coexistence 
curve  (0  = ±1), 

Ap  = ±*pc[(l -&2r>Ar/rcje,  (2) 

and  the  relation  at  constant  r, 

A p(T  = Tc)  = igpcb~'l&T(p  =pe)rr  ] \ (3) 

which  connects  the  densities  at  the  critical  tem- 
perature with  the  temperature  at  the  critical  den- 
sity for  the  same  r. 

The  dimensionless  constants  may  be  adjusted  so 
that  (i)  the  6 = ±1  curves  [Eq.  (2)]  approximate  the 
two  branches  of  the  experimental  coexistence 
curve  and  (ii)  the  T -Tc  and  p =pc  intercepts  of 
contours  of  constant  r [Eq.  (3)]  occur  at  approxi- 

AT/TC 


FIG.  2.  Contours  of  constant  susceptibility  in  C02. 
Data  are  from  Ref.  10.  The  solid  lines  are  empirically 
determined  fits  to  the  data  (Refs.  11,  19).  Note  the  ex- 
panded scales  used  in  Fig.  2(b). 
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FIG.  3.  (a)  Contours  of  constant  6 (solid  lines)  and 
constant  r (dashed  line)  in  the  T-p  plane.  The  critical 
point  is  at  r-Q.  (b)  The  form  of  these  contours  given  by 
Eq.  (1)  when  (3  = 0.35  and  g = b2  = 1.3. 

mately  the  same  locations  as  the  corresponding 
intercepts  of  constant  susceptibility  (or  constant 
correlation  length).  An  analysis  of  the  coexistence 
curve  of  C02  (Ref.  12)  and  of  contours  of  constant 
susceptibility10  " ,24  indicates  that  requirements 
(i)  and  (ii)  lead  to  J ^ 0.347  with  rc  = 31.00°C, 
g(b2  - 1 )'fl  1 .96,  g/b^\. 12,  where  in  each  case 
the  last  significant  figure  is  uncertain  by  - ±3. 

(The  critical  density  is  about  236.8  ±0.1  amagat 
= 0.468  g/cm1.)  This  results  in  the  following  set 
of  values  for  3,  g,  62: 


3 = 0.35, 

(4a) 

*=62=1.3. 

(4b) 

The  constants  (4)  have  been  rounded  off  to  two 
significant  figures  for  use  as  a first  approximation 
consistent  with  and  within  the  error  bars  of  the 
measurements  of  (2)  and  (3). 

The  prediction  (3)  for  this  choice  of  constants  is 


compared  with  susceptibilities  measured  over  a 
wide  range  of  temperature  and  pressure  in  Table 
I.  Agreement  over -all  tends  to  be  within  '3%  when 
experimental  intercepts  at  Ap£0  are  averaged, 
with  somewhat  better  agreement  for  the  Ap  > 0 
branch  than  for  the  A p<0  branch. 

The  discrepancy  between  real  and  calculated 
quantities  in  Table  I results  partly  from  errors  in 
choosing  the  constants  (4)  and  partly  from  the  fact 
thatcontours  of  constant*' in  the  linear  model  (see  Fig. 
3)  and  contours  of  constant  susceptibility  that  have 
been  constructed  from  the  experimental  data11’19,24,25 
(see  Fig.  2)  do  not  exactly  coincide  at  the  critical 
isotherm,  0 = 1/6,  as  well  as  at  the  critical  iso- 
chore, 6 = 0.  However,  the  major  discrepancy 
between  the  real  and  calculated  values  in  Table  I 
results  from  the  fact  that  the  fluid  is  inherently 
asymmetric  about  the  critical  density,17’ 19  with  the 
asymmetry  increasing  as  1 AT|  increases  (see  Fig. 
2),  while  the  linear  model  is  “everywhere”  sym- 
metric. A quantitative  approach  to  the  representa- 
tion of  this  asymmetry  is  presented  else- 
where.11’ls’ 25  For  the  purposes  of  this  paper  it 
is  sufficient  to  point  out  that  as  the  critical  point 
is  approached  the  asymmetry  becomes  less  and 
less  evident,  so  that,  within  some  empirical  bound- 
ary, the  fluid  data  may  be  represented  quite  ac- 
curately by  a totally  symmetric  equation  of  state. 
(For  C02  this  boundary  is  roughly  a contour  of 
constant  r with19  r~ 0.01.) 

Ill  COMPRESSIBILITY  AS  A FUNCTION  OF 
TEMPERATURE 

The  exponent  > can  be  estimated  independently 
of  g,  b2,  and  /I  by  considering  the  dependence  of 
the  susceptibility  on  temperature  along  the  critical 
isochore  (p -pc , 0=0).  Measurements  made  over 
a wide  range  of  AT>0  have  given  results1’’"  which 

TABLE  I.  Maximum  temperature.  Tt  —Tc.  with  7C 

31.00  C,  and  intercepts  on  the  critical  isotherm  of  con- 
tours of  constant  susceptibility  (Fig.  2):  (a)  experimen- 
tal results.  Ref.  10:  (b)  intercepts  predicted  by  Eqs.  (3) 
and  (4).  The  experimental  values  for  Ap  are  uncertain 
to  ~ 1%  of  Ap. 


(a)  Experiment 
X T,-Tc 

(102  amagat2  atm'1)  CC) 

Ap  (T  Tc  ) 
(amagat) 

(b)Theory 
A p(T  Tc\ 
(amagat  i 

100 

3.83 

-56,  +58 

± 58.5 

50 

6.67 

-67,  +73 

* 71.0 

25 

11.8 

-83.  +88 

± 86.6 

14 

19.1 

-97.  + 106 

i 102.4 

10 

26.1 

-107.  + 120 

± 114.3 

5.4 

44.1 

-126.  ♦ 144 

± 137.3 

3.3 

68.6 

-142.  + 169 

± 160.3 

2.0 

109 

-161.  ♦ 191 

± 187.8 
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are  comprised  in  the  statement 

K'Tl  ^16(0.01)1•2','F<.(AT/T<.)^  (5) 

with13  y^l.17  or14  y-  1.21 9.  The  different  values 
of  y reflect  differences  of  interpretation  of  light  - 
scattering  data,  experimental  uncertainties,  and 
a difference  in  weights  attached  to  data  obtained  at 
the  extremes  of  the  temperature  range.  In  par- 
ticular, the  lower  exponent  gives  better  agreement 
with  data  obtained  far  from  Tc  (Table  I)  and  the 
higher  exponent  better  agreement  close  to  Tc. 

(The  light -scattering  data  of  Ref.  13  are  compatible 
with  a value  of  the  compressibility  exponent  that 
varies  from  y2  1.20  for  small  AT  to  ys  1.17  for 
large  AT.) 

The  approximations 


y = l.2, 

(6a) 

n = 21 

(6b) 

are  consistent  with  the  choice  (4)  and  the  relations 
(A7),  (El),  (B6),  and  (5).  (Equations  whose  num- 
bers are  preceded  by  a letter  are  found  in  the 
Appendix  of  the  same  letter.)  More  generally, 
for  different  choices  of  y,  the  constant  a is  con- 
strained by  (B6)  and  (5)  to  have  the  value 

a = 16g(0.01)1,2">.  (7) 

In  the  model  critical  state  (Bl),  the  susceptibility 
(B6)  is  smaller  than  (5)  for  AT<0  and  d = ±l  by  a 
factor  xJx.  = 2(62  - l)l_r|l  -*2(1  -2(3)]‘‘  that  be- 
comes =4.17  for  the  choices  (4),  (6a).  The  experi- 
mentally determined  factor  is13  x+/x.~4.0,  with 
an  uncertainty  ~0.3.  The  experimental  uncertainty 
of  (5)  at  AT/TC  = 0.01,  and  hence  the  uncertainty  of 
the  numerical  coefficient  in  (7),  is  approximately 
5%. 

A choice  y = 1.2  and  a =20,  to  be  employed  in 
some  of  the  discussion  to  follow,  is  compatible 
with  (7)  within  its  error  bar.  This  choice  of  con- 
stants provides  a somewhat  better  average  fit  to 
(5)  than  does  (6)  in  the  case  that  the  empirical 
relation  (5)  actually  should  be  used  with  y varying 
from  y > 1.2  for  small  | AT  | to  y <1.2  for  large 
| aT|.  The  choice  y = 1.20  is  then  an  average  value 
ar.d  the  choice  <t=20  is  a best-fit  value  for  use 
when  the  value  y = 1.20  is  used  over  the  whole  tem- 
perature range. 

It  may  be  noted  that  the  choice  y = 1.2,  a =20, 
and  the  choice  of  constants  given  by  (4)  agrees  to 
within  6' t with  the  corresponding  parameters 
(d  = 0.3486,  ft2  - 1 .35,  g = 1.37,  a-19.6,  y = 1.178) 
found  by  Murphy,  Senger3,  and  Sengers"1  in  a fit 
to  Michels's  pressure,  density,  and  temperature 


data10  on  C02  using  the  “restricted”  linear  model. 
There  is  substantial  disagreement,  however,  with 
their  unrestricted  linear-model  results  (d  =0.3486. 
62=1.70,  g = 1.75,  a = 24.48,  y = 1 .175)  which  they 
found  to  give  better  agreement  with  the  same  pres- 
sure data.  It  is  to  be  noted  that  their  result  b2 
= 1.70  for  this  latter  mode)  is  incompatible  with  an 
interpretation  that  contours  of  constant  r should 
closely  approximate  contours  of  constant  sus- 
ceptibility (or  constant  correlation  length).  For 
this  reason,  and  others  mentioned  in  Appendix  BD 
and  in  Sec.  V,  we  prefer  to  employ  an  unrestricted 
linear  model  with  a smaller  value  of  b2.  The 
“restricted”  model  parameters  (3=0.347,  h2=1.44, 
^ = 1.49,  a = 28.37,  y = 1.24  obtained  by  Hohenberg 
and  Barmatz15  differ  more  from  our  own  than  the 
more  recent  restricted  model  set  of  parameters 
found  by  Murphy  et  al.'e 

IV.  SUSCEPTIBILITY  AND  CHEMICAL  POTENTIAL  AS 
FUNCTIONS  OF  DENSITY  AND  TEMPERATURE 

A detailed  comparison  of  the  predictions  (Blc) 
and  (B6)  with  the  choice  of  constants  ,3  = 0.35, 

# = 62=1.3,  o=20,  y = 1.20  and  experimental  data 
derived  from  bulk  measurements  of  pressure, 
volume,  and  temperature10  is  presented  in  Figs. 

4 and  5. 

Agreement  with  the  average  observed  chemical 
potential17  (Fig.  4)  is  to  within  =3%.  (Individual 
data  points  have  an  rms  deviation  .)  The  data 
correspond  to  values  of  r in  the  range  0.0007 
<r<0. 03,  for  which  there  is  almost  perfect  anti- 
symmetry about  Ap  = 0 (e  = 0).  Scatter  in  the  data 
is  greatest  for  r<  0.003,  wiic.-y  the  experimental 
uncertainties  in  AT,  Ap,  Ap  are  greatest. 

The  inverse  susceptibilities  plotted  in  Fig.  5 
encompass  a wider  domain  of  Ap  and  AT  (corre- 


FIG.  4.  Dependence  of  the  chemical  potential  on  8 (den- 
sity! at  constant  r.  Data  are  from  Ref.  17,  plotted  for 
Tc  = 31.00 °C.  The  solid  curve  corresponds  to  Eq.  (Blc). 
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FIG.  5.  Susceptibility  as  a function  of  6 and  r.  Data 
are  from  Ref.  10,  plotted  for  7>  31.00  “C.  The  solid 
line  is  given  by  Eq.  (B6). 


the  constraint  (7)  is  satisfied,  other  choices  of 
y in  the  range  y = 1.20  ±0.03  yield  results  similar 
to  those  shown  in  Figs.  4 and  5. 


V SPECIFIC  HEAT 


The  experimental  dependence  on  temperature  of 
the  specific  heat  of  C02  is20 


(8) 


where  R is  the  gas  constant  and  A*  =5.583,  A. 

= 10.473,  B,  =-3.457,  B.  = -0.024,  a*  = a,  =0.125. 
(The  subscripts  plus  and  minus  refer  to  tempera- 
tures, respectively,  above  and  below  7’c.) 

The  experimentally  determined  ratio 

A. /A*  =1.876  (9) 


may  be  compared  with  that  predicted  for  the  model 
critical  state  [Eqs.  (B5)  for  0=0  and  (B9)j ; 


sponding  to  r«0.17)  and  show  convergence  toward 
linear -model  behavior  [Eq.  (B6)l,  which  is  rep- 
resented by  the  solid  line,  as  r decreases.  For 
the  restricted  range  r<0.03  used  in  Fig.  4,  the 
susceptibility  remains  constant  to  approximately 
5%  along  most  of  a contour  of  constant  r.  The 
rather  good  agreement  between  the  predictions  of 
the  simple  (linear)  parametric  model  and  the  sus- 
ceptibility data  is  not  unexpected,  because  a care- 
ful study  of  these  data  by  themselves"  has  led  us  to 
conclude  that  a parametric  equation  is  a “natural” 
equation  for  the  representation  of  such  data  (see 
Appendix  C). 

Appreciable  deviations  from  constancy  and  from 
(B6)  are  confined  primarily  to  larger  values  of 
r and  the  regions  |0|>O.8,  for  which  |Ap|  is  large. 
The  largest  deviations,  for  0<-O.8,  reflect  the 
presence  of  the  absolute  boundary  at  Aplm„  = -pc 
(rarified  gas)  at  which  — The  behavior  of  the 
susceptibility  near  this  boundary  is  expressed 
most  naturally  in  terms  of  p,  T rather  than  r,  6 
coordinates;  the  divergence  of  X'1  in  this  region 
of  dilute  gas  corresponds  to  a component  of  (x"')& 
in  (El)  proportional  to  p"1.  The  deviations  from 
the  prediction  (B6)  for  large  r and  0>+O.8  may 
correspond  to  a background  (x*')6  that  continues  to 
be  dependent  on  p when  p is  large. 

Additional  deviations  for  small  r and  |S|  = 1 can 
be  expected  if  the  anomalous  behavior  near  | 0| 

= |0m|  and  its  extreme  sensitivity  to  the  choice  of 
value  for  ft2  (cf.  the  discussion  in  Appendix  B D) 
are  indicative  of  a defect  of  the  model  equation  of 
state  or  choice  of  ft2  that  we  have  used.  Any  such 
additional  deviations  for  small  r seem  at  present, 
however,  to  fall  within  the  limits  of  accuracy  of 
existing  data. 

For  the  same  choices  of  h2,  (3,  and  provided 


^4-  _ /n  + J?  + 1 4 

A.  ■ (h2-l)>‘26/„  ' 


(10) 


For  62  = 1.30,  (1=0.35,  the  theoretical  ratio  is 


A_/A+  = 2.07  for  y = 1.20,  (11a) 

A_/At  = 1.76  for  y = 1 .22.  (lib) 


The  predictions  (11)  are  sensitive  to  (1  (rather  in- 
sensitive to  I)2).  The  choice 

y = 1 .217,  (3=0.347,  g = h2  = 1.30,  (12) 


which  gives  A ./A,.  =1.86,  is  compatible  with  the 
experimental  ratio  (9)  and  also  with  the  shape  of 
the  coexistence  curve  (Sec.  II).  Significantly 
poorer  agreement  with  the  experimental  ratio  of 
specific  heats  is  obtained  using  either  the  param- 
eters of  the  restricted  linear  model  or  those  of  the 
unrestricted  linear  model  employed  by  Murphy 
el  al.'6  Using  the  parameters  quoted  by  Murphy 
cl  al.,  the  theoretical  ratio  (10)  becomes  A. /A* 

= 2.54  (restricted  model)  or  2.98  (unrestricted 
model).  The  theoretical  value  for  the  ratio  ob- 
tained by  using  the  restricted-linear-model  pa- 
rameters of  Hohenberg  and  Barmatz'5  is  A.,/A» 

= 1.53.  Use  of  our  parameter  values  ) =1.217  and 
,3=0.347  in  the  Ho-Litster"  (cubic)  model  (cf.  end 
of  Appendix  C)  leads  to  values  of  b2  =fe2,,  =3/(3  - 2,1) 
and  of  the  specific -heat  ratio25  almost  identical  to 
the  values  (12)  and  (9),  as  found  above  using  our 
version  of  the  unrestricted  linear  model. 

Above  Tc  along  the  critical  isochore  (9  = 0)  the 
specific  heat  (B5)  is 


(fP).  = 


-t>HPcc0T 

T\ 


(xV"- 


(13) 


This  expression  cannot  exactly  reproduce  the 
temperature  dependence  stated  in  (8)  when  2 -> 
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-2*1  = <**0.125.  However,  as  discussed  in  Ap- 
pendix D,  for  values  of  <*  not  too  far  different 
from  0.125  there  is  good  numerical  agreement 
for  suitable  choices  of  a and  g.  For  example,  for 
£ = 1.30  and  *1  = 0.35,  agreement  of  (13)  with  the 
temperature -dependent  part  of  (8)  above  Tc  is  to 
within  ~1%  over  most  of  the  range  of  temperatures 
for  which  reliable  data  are  available  for  the  choices 


y = 1.20, 

<*=0.10, 

a = 23.4, 

(14a) 

y = 1.22, 

<*  =0.08, 

a =23.9. 

(14b) 

Agreement  with  Eqs.  (B9)  and  (B8)  below  Tc  to 
~1%  is  obtained  for  the  same  choices  of  y and  <* 
for  somewhat  different  values  of  a.  Table  II  sum- 
marizes values  of  a,  and  added  background  Bt 
that,  used  with  Eqs.  (13)  and  (B9),  fit  the  specific - 
heat  data  (8)  above  (+)  and  below  (-)  Tc  to  '1%  for 
the  choices  of  exponents  discussed  above.  In  the 
same  table  are  listed  the  values  of  a required  to 
fit  the  susceptibility  l(B6a),  (7)]  for  the  same 
choices  of  exponents. 

When  the  scaling-law  requir  ner  that  B+  and 
B_  be  equal17  is  applied  for  j = b2=  1.3, 

y=1.20,  a =20,  and  B.  = B_  = -7  there  is  over-all 
agreement  with  both  branches  of  the  specific  heat 
and  with  the  compressibility  (7)  to  better  than  10%. 
Over-all  agreement  improves  for  each  of  the  fol- 
lowing modifications:  (a)  slightly  smaller  *1,  (b) 
somewhat  larger  y compatible  with  (7).  As  an 
example,  the  choices  g=b2  =1.30,  |3  = 0.347, 
y = 1.217,  n=22  [Eqs.  (7),  (12)j  used  with  equal 
background  contributions  B . = B_  = -8.9  lead  to 
agreement  with  all  of  the  specific -heat  data20  and 
the  compressibility  to  =3%. 

Agreement  with  the  specific  heat  in  the  range 
10‘4<|(r  - Tc)/Tc]<  10'2  measured  both  above  and 
below  Tc  can  be  improved  to  ~1^  without  noticeably 
changing  the  agreement  with  the  compressibility 
with  further  small  modifications  of  the  values  of 
the  parameters  used  in  the  equations.  Examples 
of  such  modifications  that  lead  to  improved  agree- 
ment with  (8)  for  the  same  choice  of  background 
(B . - -8.9)  include  use  of  *1  = 0.346,  a =0.087, 

62  = 1.31;  use  of  *1  = 0.344,  <*=0.089,  and  either 
g - 1 .28  or  b2  = 1 ,31;  or  use  of  *1  = 0.342,  <*=0.090, 


# = 1.28,  62=1.31,  « = 22.2.  There  remains  a small 
possibly  significant  discrepancy  with  (8),  especial- 
ly for  small  |A 7"|.  The  discrepancy  grows  to  >2% 
for  (T -Tc)/Tc=± 4X10'5,  the  smallest  A T studied 
experimentally.  This  discrepancy  cannot  be  re- 
moved by  making  additional  small  changes  in  the 
values  of  the  parameters.  Further  consideration 
of  this  point,  which  will  not  be  undertaken  here, 
as  well  as  efforts  to  improve  significantly  on  the 
values  of  the  parameters  cited  in  the  previous 
paragraph,  should  take  into  consideration  effects 
of  compression  of  fluids  under  their  own  weight 
and  the  lengthening  of  thermal  equilibration  times 
that  occurs  close  to  the  critical  point.  Considera- 
tion might  also  have  to  be  given  to  the  possible 
effect  on  the  numerical  values  of  these  param- 
eters caused  by  contamination  or  whatever  other 
factor  accounts  for  the  difference  between  the 
values  of  Tc  found  in  most  of  the  experiments  dis- 
cussed (31.0°C)  and  that  obtained  in  the  calorimet- 
ric measurements  (30.8  C).20 

This  completes  our  analysis  of  experimental 
data  to  obtain  a set  of  constants  to  use  in  the  un- 
restricted-linear-model equations  of  the  critical 
state . The  more  refined  of  the  two  sets  of  parameters 
with  g = 6s  = 1.30  discussed  above  (that  accurate  to 
-3%)  is  listed  in  Table  III  together  with  values  for 
these  parameters  obtained  in  other  recent  investi- 
gations. 

The  constants  used  in  the  model  equation  of  state 
have  been  determined  above  from  experiments 
performed  sufficiently  far  from  the  critical  point 
that  the  compression  of  a fluid  under  its  own 
weight  is  of  secondary  importance.  In  the  follow- 
ing sections  the  model  thus  specified  will  be  used 
to  interpret  experiments  performed  to  study  ef- 
fects of  this  compression.  The  results  of  these 
experiments  will  be  examined  to  see  if  they  verify 
the  theory  for  regions  of  thermodynamic  space 
that  are  too  close  to  the  critical  point  to  be 
studied  effectively  by  the  types  of  investigations 
considered  to  this  point. 

VI.  DENSITY  AS  A FUNCTION  OF  HEIGHT 

According  to  (F2),  (E2),  and  (Bl),  the  density 
of  a linear  model  fluid  at  constant  temperature 


TABLE  II.  Values  of  a and  B required  to  fit  the  measured  compressibility  l Eq.  (7)1  and 
specific  heat  1 Eqs.  (8),  (B9>.  and  (13)1  forg  t>2  1 .30  and  different  choices  of  critical  expon- 

ents. Subscripts  plus  and  minus  refer  to  temperatures,  respectively,  above  and  below  Tc. 


Exponents 

Compressibility 

Specific  heat 

o 

0 

y 

a t 

a _ 

B ¥ 

B _ 

0.10 

0.35 

1.20 

21 

23.4 

21.2 

-7.1 

-7.0 

0.08 

0.35 

1.22 

23 

23.9 

25.5 

-11. S 

-15.6 

0.089 

0.347 

1.217 

22.7 

22.9 

23.2 

-9.4 

-11.4 
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TABLE  lit.  Parameters  for  use  with  the  “linear-mod- 
el" representation  of  the  critical  region  of  C02.  Columns 
marked  with  an  asterisk  correspond  to  the  “restricted” 
model,  for  which  b2  is  constrained  to  satisfy  Eq.  (BIO). 


Source 

Hohenberg 
and  Barmatz** 

Murphy 

etal** 

Murphy 
et  al  . b 

Present 

investigation 

Q 

0.065 

0.13 

0.13 

0.089 

p 

0.347 

0.3486 

0.3486 

0.347 

Y 

1.241 

1.178 

1.175 

1.217 

a 

28.4 

19.6 

24.5 

22.0 

b2 

1.44 

1.35 

1.70 

1.30 

g 

1.49 

1.37 

1.75 

1.30 

a Reference  15. 


Reference  16. 


varies  with  height  in  accordance  with  the  para- 
metric relation 

A z=azytR(e-0'),  (15) 

with 

-z0=Pc(mpcM0)'1.  (16) 

Here  Az  is  the  distance  above  that  elevation  z(pc) 
in  the  fluid  at  which  p-pc.  It  follows  from  (Bl) 
and  (15)  that  for  T>  Tc  and  small  Ap 


Az  = 


= znaAp  / AT V 

gPc  v Te  I ’ 


and  for  T = TC,  Ap*0, 

z.u&s*3(6 2 - l)Ap 

Ac  = r — 

8 Pc 


Ap 

Pc 


1 


(17a) 


(17b) 


where  6 = (y/ji)  + 1 . 

The  predictions  (17)  are  in  qualitative  agree- 
ment with  density  profiles  measured  in  C02  by 
Schmidt  and  Straub.21'22  Quantitatively,  however, 
the  value  of  a required  to  fit  Schmidt  and  Straub’s 
data  is  much  smaller  (by  about  a factor  of  3)  than 
the  value  (6). 

The  prediction  (15)  Imore  generally,  the  exact 
proportionality  (FI)  between  increments  of  eleva- 
tion and  chemical  potential]  depends  crucially  on 
the  absence  of  thermal  gradients.  When  thermal 
gradients  are  present,  it  is  necessary  to  use  (F5) 
in  place  of  (FI).  In  C02,  for  which  at  the  critical 
point17-22  (8/yar)p  = 6.99/Jc  T~'  = 1.7  xlO6  erg/cm3 
deg,  a gradient  7,=  <-0.6  mdeg/cm  (top  of  the  cell 
warmer  than  the  bottom)  could  account  for  the  dis- 
crepancy between  the  observed  density  profile  and 
that  predicted  by  (15).  A thermal  gradient  of  this 
magnitude  appears  to  have  been  within  the  limits 
of  accuracy  of  the  temperature  measurements  in 
the  experiments.21 

VII.  SUSCEPTIBILITY  AND  SCATTERING  PROFILES 

Expressions  (Bib),  (B6),  (15)  can  be  solved  to 
find  the  dependence  of  the  susceptibility  of  the 


ideal  critical  fluid  on  height  at  constant  tempera- 
ture. For  A7>0,  the  prediction  to  first  approxi- 
mation in  (A zY  is  (ignoring  the  small  dependence 
of  x on  S2  at  constant  r) 


X«  = 


8Pc 

aPr 


/A7\-’r 

>i>2 

/Ary2(),+8,/A2Y7 

\TJ  L1' 

a2 

[tJ  \ zj  } 

(18) 


Since  the  intensity  of  light  scattered  from  each 
element  of  a fluid  is  proportional  to  the  suscepti- 
bility,1"3 this  leads  to  the  prediction  of  scattered 
light  profiles  of  approximately  the  form  (18). 

The  prediction  (18)  for  the  choice  of  constants 
(4),  (6)  is  compared  with  the  authors'  light -scat- 
tering measurements23  in  Fig,  6.  The  solid  curves 
are  experimental.  Each  dashed  curve  corresponds 
to  Eq.  (18)  multiplied  by  a constant  scale  factor  to 
bring  the  top  of  the  curve  into  coincidence  with 
the  peak  scattered  intensify  at  the  same  tempera- 
ture. A smaller  scale  factor  is  required  for  the 
upper  curve  (smaller  AT)  than  for  the  lower  curve. 

The  need  for  different  scale  factors  and  other 
discrepancies  between  the  susceptibilities  given  by 
(18)  and  the  measured  scattered  intensities  are 
qualitatively  those  to  be  expected  to  result  from 
the  partial  extinction  of  incident  and  scattered 
light  which  occurs  in  a finite  thickness  of  fluid 
near  the  critical  point.  The  extinction  causes  a 
reduction  of  observed  scattered  intensity.  The 
fraction  of  the  incident  and  scattered  light  re- 


EIG.  6.  Scattered  intensity  l and  susceptibility  % as 
functions  of  height.  The  solid  curves  give  intensities 
observed  at  13.5°  scattering  angle;  each  dashed  curve 
corresponds  to  Eq.  (IS),  multiplied  by  an  arbitrary  con- 
stant. See  text. 
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moved  by  extinction  increases  as  the  critical  point 
is  approached  more  closely.  Consequently,  within 
an  isothermal  sample  the  extinction  increases  to  a 
maximum  value  that  is  characteristic  of  the  par- 
ticular temperature  as  the  “critical  height”  (2  = 2e, 
p ~pc)  is  approached  either  from  above  or  below. 
The  effect  of  the  height -dependent  extinction  on  the 
scattering  profiles  observed  in  a sample  of  fluid 
of  finite  thickness  is  therefore  a decrease  in  the 
® amplitude  of  the  top  of  each  profile  relative  to  the 
amplitude  of  the  wings,  combined  with  an  over-all 
reduction  of  height  which  is  dependent  on  tempera- 
ture. These  features  are  apparent  in  Fig.  6 which, 
in  effect,  compares  profiles  measured  for  a sam- 
ple of  finite  thickness  (solid  lines)  to  profiles  that 
would  be  predicted  for  an  ideally  thin  sample  of 
the  same  fluid  (dashed  lines)  if  the  amplitudes  of 
the  latter  profiles  are  “corrected”  by  temperature - 
dependent  factors  to  bring  the  peak  intensities  into 
coincidence  with  those  of  the  observed  profiles. 

The  need  for  the  temperature -dependent  correc- 
tion factors  is  also  apparent  in  the  figure,  which 
shows  that  the  ratio  of  observed  peak  intensities 
(solid  curves)  differs  by  a factor  less  than  2 land 
considerably  less  than  predicted  by  Eq.  (18)]  for  a 
larger -than-2  ratio  in  the  temperature  distances 
from  the  critical  point. 

Quantitatively,  the  magnitude  of  the  corrections 
for  extinction,  as  also  those  for  (a)  angular  ani- 
sotropy of  the  scattering  and  (b)  multiple  scatter- 
ing, depends  nonlinearly  on  the  magnitude  of  the 
susceptibility  and  on  the  correlation  range.  A 
fuller  discussion  of  these  corrections  may  be  found 
in  Ref.  13,  where  an  attempt  was  made  to  calculate 
them  on  the  basis  of  estimates  for  the  magnitude 
of  the  susceptibility  and  the  correlation  length 
close  to  the  critical  point. 

VIII  CONTOURS  OF  CONSTANT  WAVE  SCATTERING 

Equations  (Bib)  and  (15)  define  temperature  and 
height  coordinates  of  contours  of  constant  sus- 
ceptibility (B6).  Typical  contours  are  shown  in 
Figs.  1(d)  and  7(a).  These  contours  are  approxi- 
mated rather  closely  by  contours  of  constant  r 
I Fig.  7(b)  1 . They  can  be  expected  also  to  approxi- 
mate contours  of  constant  correlation  length  (. 
Wave  scattering  from  each  element  of  the  fluid  and 
corrections  for  finite  sample  thickness  depend 
principally  on  the  range  of  correlation  and  the  sus- 
ceptibility at  each  element.  Therefore  it  can  be 
expected  that  to  a good  approximation  contours  of 
constant  scattering  from  elements  of  the  fluid  are 
simultaneously  contours  of  constant  corrections 
to  scattering,  and  hence  should  coincide  with  con- 
tours of  constant  observed  scattering,  obtained 
from  samples  of  finite  thickness.  Conversely,  it 


can  be  predicted  that  contours  of  constant  observed 
scattering  should  be  contours  of  constant  correc- 
tions to  scattering,  and  hence  contours  of  constant 
susceptibility  and  constant  correlation  range. 

This  prediction  is  compared  with  recent  experi- 
mental results2’  in  Figs.  8(a)  and  8(b).  The  con- 
tours shown  are  the  contours  of  constant  )'  pre- 
dicted by  Eqs.  (Bib)  and  (15)  for  the  choice  of 
constants  (4),  y = 1.20,  a=20,  and  7'c  = 31.0027°C. 
The  data  have  been  taken  from  individual  isother- 
mal profiles  and  refer  (Fig.  8(a)]  to  the  inten- 
sity” / of  light  scattered  through  an  angle  of  13.5°, 
and  [Fig.  8(b)]  to  the  inverse  of  the  Rayleigh  line- 
width26  (fluctuation  decay  rate)  of  this  scattered 
light.  It  will  be  noticed2’  that  to  within  the  preci- 
sion of  these  measurements,  the  intensity  contours 
and  contours  of  constant  fluctuation  decay  times 
are  of  the  same  form.  This  confirms  by  light- 
scattering  measurements  a similarity  of  scaling - 
law  behavior  for  equilibrium  and  nonequilibrium 
properties  demonstrated  by  Sengers  and  Keyes  in 
their  analysis  of  thermal -conductivity  data  for 
CO,.27 


FIG.  7.  (at  Contour  of  constant  susceptibility  (r  = const! 
as  a function  of  temperature  ST  and  height  Az.  (b)  Con- 
tours of  constant  r (dashed  lines!  and  constant  0 (solid 
lines!  for  Tc  = 304  *K,  ft2  = 1.3,  a = 21,  *0  = 1.61*  10s  cm, 
and  y + ft=1.55  lEqs.  (Bib!  and  (15)1. 
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FIG.  8.  Contours  of  constant  wave  scattering:  (a)  in- 
tensity; (b)  inverse  Rayleigh  linewidth.  Solid  curves  are 
contours  of  constant  r,  plotted  for  Tc  =31.0027 °C.  Val- 
ues of  r have  been  assigned  arbitrarily  in  (a)  and  have 
been  set  equal  to  (t/t0)'1'1'  in  (b) , with  r0  = 1.55  psec, 
v = 0.63. 


I 


: 


The  contours  of  constant  r shown  in  Figs.  7(b), 
8(a),  and  8(b)  have  an  aspect  ratio 


Az,„lx 

at;,,,. 


-arr*R~' 


z0  Tc 


(19) 


Here  Az,,,,,  is  the  maximum  value  of  a -zc  and 
AT„ax  is  the  maximum  value  of  T - Tc  for  each 
contour.  [Equation  (19)  follows  immediately  from 
Eqs.  (Bib)  and  (15).]  The  choice  of  parameters 
(7)  and  (12)  leads  to  slightly  smaller  aspect  ratios 
than  those  shown  in  Figs.  8(a)  and  8(b).  Suscepti- 
bility contours,  derived  from  Eq.  (B6),  have  a 
slightly  different  shape  (longer  and  narrower)  and 
slightly  smaller  aspect  ratio  than  contours  of  con- 
stant r.  The  combined  effect  of  the  two  corrections 
[use  of  Eq.  (B6)  and  parameters  (7)  and  (12)] 
amounts  to  — 109t  deviations  from  the  contours  ex- 
hibited in  Figs.  8.  Existing  data,  shown  in  Fig.  8, 
confirm  the  general  correctness  of  the  theory  and 
of  the  estimates  of  parameters  arrived  at  from  a 
consideration  of  other  data,  mostly  obtained  far- 
ther from  the  critical  point,  but  are  insufficiently 
precise  to  lead  to  a definitive  test  or  refinements 
of  those  results.  An  order-of-magnitude  improve- 
ment of  experiments  of  this  type  appears  readily 
attainable,  however,  and  could  lead  to  significant 
new  tests  of  the  behavior  of  a fluid  close  to  its 
critical  point.  Experiments  directed  to  this  end 
are  currently  in  progress,  and  will  be  reported 
elsewhere. 


IX.  CONCLUSIONS 

We  have  argued  that  in  carbon  dioxide  an  ex- 
pression for  the  free  energy  of  the  form  (E3), 

(B3)  can  be  found  which  predicts  correctly  the 
magnitude  of  the  critical  components  of  all  bulk 
thermodynamic  properties  (compressibility,  spe- 
cific heat,  chemical  potential)  and  which  is  com- 
patible with  existing  light-scattering  measure- 
ments. In  particular,  no  evidence  is  found  to  with- 
in the  precision  of  available  experimental  data  and 
the  accuracy  of  this  analysis  for  any  departure 
from  the  fundamentally  nonanalytic  behavior  of  the 
thermodynamic  functions  close  to  the  critical  point 
and  the  simplicities  of  their  behavior  that  are 
expressed  in  the  recently  discussed  general  laws 
of  scaling,  homogeneity,  and  universality .28",° 

Quantitative  comparisons  of  the  simple  formu- 
lation (Bl),  (B3)  that  incorporates  these  laws  with 
experimental  data  have  shown  that  for  the  (first- 
approximation)  choice  of  the  five  dimensionless 
constants  that  appear  in  the  model  expression  (B3) 
for  the  free  energy,  namely.  62=g  = 1.3,  a=21. 
(3=0.35,  y = 1.20,  and  for  a constant  background 
contribution  to  the  specific  heat,  agreement  with 
most  of  the  measurements  in  the  range  3xi0"5 
<r<3><10"2  is  within  ~10%.  Improved  agreement 
(~5%)  is  obtained  for  the  choice  J = 0.347.  y = 1.217, 
a =22,  and  a constant  background  for  the  specific 
heat.  Agreement  is  to  better  than  5%  with  both 
branches  of  the  specific -heat  curve  and  with  some 
of  the  most  accurate  PVT  data  close  to  the  critical 
point.  Disagreement  tends  to  be  greatest  at  the 
largest  temperature  and  density  distances  from 
the  critical  point,  where  perfect  scaling  sym- 
metries begin  to  break  down.  Further  improve- 
ment can  be  obtained  by  employing  background 
contributions  that  change  gradually  with  tempera- 
ture and  density  or  by  making  small  changes  in  the 
form  of  the  model  equations  of  state,11  •18,lS1’25 
possibilities  that  have  not  been  analyzed  in  detail 
here.  For  the  model  (Bl),  (B3)  and  either  of  the 
above  choices  of  constants,  the  susceptibility  and 
the  correlation  length  vary  by  only  a few  percent 
over  the  full  range  of  -1  =s0s>+l  for  constant  r; 
hence  this  the  “distance”  from  the  critical 
point,  is  in  the  model  used  here  a measure  to  a 
first  approximation  of  the  mean  range  of  local 
correlations  and,  via  the  equivalence  of  contours 
shown  in  Figs.  8(a)  and  8(b),  also  a measure  to  a 
first  approximation  of  the  mean  time  for  growth 
and  decay  of  fluctuations  of  localized  ordering. 

The  numerical  values  of  the  exponents  considered 
here  differ  significantly  from  the  values  ,1=0.3125, 
y=1.25  predicted  by  the  Ising  model  of  the  critical 
region.  It  may  be  noted  that  a value  of  y<  1.20 
is  found  when  experimental  data  obtained  close  to 
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the  critical  point  are  ignored  or  assigned  a re- 
duced weight.  A possible  increase  of  y to  y > 1 .22 
and  simultaneous  decrease  of  ji  to  appreciably 
less  than  0.35  very  close  to  the  critical  point  may 
have  been  missed  in  the  experimental  investiga- 
tions analyzed  in  this  paper;  questions  of  this 
kind  may  be  resolved  through  refinements  of  the 
technique  and  analysis  of  measurements  of  the 
type  described  by  Eqs.  (18)  and  (19)  of  this  paper. 
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APPENDIX  A:  GENER  AL  THERMODYNAMIC 
RELATIONS  AND  DEFINITIONS 

For  convenience  of  reference  we  summarize 
here  the  definitions  and  derivations  of  the  thermo- 
dynamic quantities  required  for  the  present  in- 
vestigation. 

At  thermodynamic  equilibrium,  the  energy  E of 
a fixed  quantity  of  fluid  at  temperature  T , entropy 
S,  pressure  P,  volume  V,  obeys  the  second  law 
of  thermodynamics  in  the  form 

dE=TdS -PdV.  (Al) 

As  a consequence,  its  Gibbs  free  energy, 

G = £-TS+PF=pN  (A2) 

(N  is  the  number  of  particles,  p the  chemical  po- 
tential per  particle),  is  a function  of  the  intensive 
variables  T and  P, 

dG=Ndp=-SdT  + VdP,  (A3) 

and  its  Helmholtz  free  energy  per  unit  volume. 

/=F/F=e-Ts=pp-P.  (A4) 

(e  = E/V,  s=S/V,  p-N/V),  is  a function  of  density 
and  temperature, 

dj  - p dp  ~sdT . (A5) 

Expression  (A5)  for  the  Helmholtz  energy  follows 


from  the  Gibbs -Duhem  relation  (A3),  which  re- 
quires 

cU’  = prfp  +sdT.  (A6) 


The  Helmholtz  energy  can  be  used  to  generate  all 
thermodynamic  functions  of  interest. 

First  derivatives  of  the  Helmholtz  free  energy 
yield  expressions  for  the  chemical  potential, 


(A7a) 


(A  7b) 


as  functions  of  temperature  and  density.  Higher 
derivatives  give  the  isothermal  compressibility, 


Kt 


1 (*X.\  - _L  (1 £ _ J_  ( ap) 

v\dP/T  p \ap/T  p2\dp) T’ 


(A7c) 


or  susceptibility, 

(inl- 
and the  specific  heat  at  constant  density, 


(A7d) 


(A7e) 


The  latter  two  quantities  are  of  particular  interest 
because  they  both  diverge  at  the  critical  point 
(p=pc,  T -Tc)~ the  specific  heat  weakly,  the  sus- 
ceptibility (compressibility)  strongly. 


APPENDIX  B:  THERMODYNAMICS  OF  LINEAR  MODEL 

The  general  parametric  equations  of  state  of 
the  linear  model  of  Schofield6  written  in  a form 
applicable  to  fluids  are 


Ap=gpcfa  0, 

(Bla) 

AT  = Tcr(l  -b2P), 

(Bib) 

Ap  =aPcp~c'rR*y(6  - 03), 

(Blc) 

where  pc.  Tc,  Pc  are,  respectively,  the  density, 
temperature,  and  pressure  at  the  critical  point. 
Here  Ap  =p  -pc,  AT  = T - Tc;  and  Ap  = p(p,  T) 

- p(pc,  T ) is  the  anomalous  part  of  the  chemical 
potential.  The  quantities  g,  )},  b 2,  a,  and  y are 
dimensionless  constants,  with  62>1. 

In  the  restricted  linear  model,7  b2  is  related  to 
the  exponents  ji  and  y (Appendix  BC).  In  the  un- 
restricted linear  model,  b2  is  an  adjustable  con- 
stant. The  exponents  (3,  y,  and  6 = (y/ji)  + l are 
particularly  significant  in  that  they  characterize 
the  shape  of  the  coexistence  curve  (0),  the  degree 
of  the  critical  isotherm  (6),  and  the  divergence 
of  the  isothermal  compressibility.  £r,  along  the 
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critical  isochore  (y).  Certain  characteristics  of 
the  linear  model  are  illustrated  pictorially  in 
Figs.  1 and  3.  Of  particular  interest  are  contours 
of  constant  inverse  susceptibility,  x*1  [Eq.  (A7d)|, 
which,  for  special  choices  of  62,  are  approxi- 
mately coincident  with  contours  of  constant  r in 
the  immediate  vicinity  of  the  critical  point  Isee 
Sec.  II  and  Eq.  (B6)J.  A typical  constant  x"1  con- 
tour as  seen  when  viewed  along  the  three  axes  in 
P-p-T  space  is  illustrated  in  Fig.  1.  Figure  3 
illustrates  contours  of  constant  r and  constant  8 
projected  onto  the  T -p  plane  for  a particular  choice 
of  g and  6 2 as  well  as  the  particular  values  of  8 
that  yield  the  critical  isochore  (6=0),  the  critical 
isotherm  (6  = 1/6),  and  the  coexistence  curve 
(6  = 1).  In  particular,  the  coexistence  curve  (6  = ±1) 
is  of  the  form  Ifrom  (Bla),  (Blb)J 

coexistence  boundary:  Ap  = ±#pc|(l  -b2)~'^T/Tc\6 . 


Equations  (Bla),  (Bib),  and  (B2)  are  identical  to 
Eqs.  (la),  (lb),  and  (2)  of  the  text. 

A.  Single -phase  region 

It  is  readily  verified  from  Eq.  (A 7a)  (Appendix 
A)  that  Eq.  (Blc)  corresponds  to  a contribution  to 
the  free  energy  of  the  form15 


fa  =agPerr*26f  (8), 


m = £/,,6',=/o  + /262  + /4e4.  (B3b) 

(A  complete  expression  for  the  free  energy  con- 
tains in  addition  background  terms,  discussed  in 
Appendix  E).  The  dimensionless  constants  /„  in 
(B3b)  are  functions  of  62,  0,  y (Table  IV).31  It  is 
to  be  noted  that  the  free  energy  (B3)  is  analytic 
everywhere  in  the  single -phase  region  except  at 
the  critical  point  (r  = 0),  where  there  is  an  essen- 
tial singularity. 

Other  thermodynamic  quantities  of  interest  can 


TABLE  IV.  Coefficients  for  use  in  Eqs.  (B3)-(B6). 
Here  a = 2 -y- 20  and  d2  6 2(1  — 2/3 ) . 


26*0(1—  o)(2-o) 
y — 20  -d2o 


^2  2 62o  (1  — o) 


s0=  (2  — o)/0 


’ 2 9 nh ‘ 


c.  = (l-a)st 


c2=  (y-  l)s2 


be  derived  from  (B3)  using  Eqs.  (A7)  (Appendix  A). 
The  results  are  as  follows: 

s,  = -atffJcr;1r,*2#',s(e),  (B4a) 

T~'(c„)m  = -agPcT-e2r*t2i-2c(6),  (B5a) 

(X'l).=aPcg-'p-2r',x-'(  6),  (B6a) 


s(6)  = so  + s20s,  (B4b) 

c(6)  = (1  - d*92)~l(c0  + r202),  (B5b) 

X”'(6)  = (1  - d262)-'(l  *X2',62+x;165).  (B6b) 

Here  d2  =62(1  -2 0)  and  the  coefficients  s„,  cn , 

X71  depend  on  62,  0,  y (Table  IV). 

From  (E5)  (Appendix  E),  (Blc),  and  (B3)  it  fol- 
lows that  the  pressure  contains  a nonanalytic  part 
of  the  form 

Pa  - (Pc  + Ap)  Pa  —/a 

= aPcry*6 (8  - 03)  -agPcr2*2R[f0  + (/2  - 1 )62 

+ (/,  +1)64J. 


The  single-phase  region  is  bounded  by  a two- 
branched  contour  0=-l  (gas)  and  6 = +l  (liquid). 

On  this  contour,  at  constant  temperature,  it  is 
evident  that  (p0),.„  =(pa)lKt,IKl  and  (P0)rj,  =(f>a),„„k.  • 
Hence  this  contour  is  the  boundary  of  coexisting 
gas  and  liquid  phases. 

B Two -phase  region 

Inside  this  boundary,  in  the  region  of  coexisting 
liquid  and  gas  phases,  expressions  (Bl)  apply  to 
each  phase  separately,  with  6 = * 1 for  the  liquid 
and  6 = — 1 for  the  gas.  In  this  region  the  chemical 
potential  (A2)  is  a function  only  of  temperature 
and  is  equal  to  its  value  at  the  coexistence  curve. 
In  the  linear  model.  Ap  =0  at  6 = ±1.  Hence  Ap  =0 
everywhere  in  the  two-phase  region  ( in  agreement 
with  Eq.  (E2)  of  Appendix  E if  Xj1  -0  for  this 
region],  and  /„  in  this  region  is.  from  (B3)  and 
(Bib), 

(\T  \J*2A 

(1  . h2)T'  J {fo* /,*/<)■ 


Successive  differentiations  of  (B8)  lead  to  ex- 
pressions for  the  entropy  and  specific  heat  in  the 
two-phase  region.  Of  particular  interest  is  the 
specific  heat, 


X = — 3 * ft 2 (—  1 * 2y  * 20)  x7'-62<3-2y  — 20) 


(^) 


(>  +2Mr  + 2,J-i) 


i 
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The  specific  heat  is  greater  in  the  two-phase  re- 
gion than  in  the  single -phase  region  at  the  same 
temperature  close  to  the  phase  boundary.  (By 
contrast,  the  entropy  does  not  change  discon  - 
tinuously  at  the  boundary  of  the  region  of  coexisting 
phases;  it  has  the  same  value  in  the  two-phase 
region  as  in  a single  phase  at  the  same  tempera- 
ture at  0 = ±1.) 


C.  Restricted  linear  model 


Some  special  choices  of  the  constants  62,  ,3,  y 
merit  comment.  For  p + y = i and  yb2  = |,  the  de- 
pendence on  0 of  the  susceptibility  (B6b)  dis- 
appears, as  does  that  of  the  specific  heat  (B5). 
[For  /3+y  = |,  the  isothermal  integration  of  Eq. 
(C2),  Appendix  C,  yields  Eq.  (C4)  exactly,  with 
y62  = |.8,11|  The  specific  heat  continues  to  be  in- 
dependent of  8 for  ii  + y*i  for  the  “restricted” 
linear  model,  in  which  b = 6* , where7 


(b'f 


y -2,3 6-3 

y (1  - 2,3)  ' (6-l)(l  - 2,3)  ’ 


(BIO) 


with  6 = (y/|3)  + 1 . 

The  discussion  elsewhere  in  this  paper  does  not 
assume  either  ;3  + y = § or  6=6*;  these  relation- 
ships are  satisfied  approximately  but  not  exactly. 


D.  Topography  of  parametric  models: 
metastable  states  and  spinodal 

For  T <TC,  Eqs.  (1)  lead  to  values  of  |Ap| 
that  are  never  smaller  than  those  which  corre- 
spond to  I 6|=  I 6m|  = d“' = [62(1 -2rl)|_L,2t  It  is  seen 
from  Eqs.  (B5)  and  (B6)  that  on  the  boundary  9 = 8m, 
for  particular  choices  of  62,  ;3,  and  y,  both  the  in- 
verse susceptibility  and  the  specific  heat  diverge 
and  change  sign.  Divergence  of  the  inverse  of  the 
susceptibility  (infinite  rigidity)  is  not  an  observed 
behavior  of  accessible  states  of  real  fluids.  For 
other  choices  of  62,  ,3,  and  y compatible  with  ex- 
periment, the  susceptibility  diverges  before  |0| 
reaches  |0„|,  at  a boundary  |0|  = |0„|,  with  |0„|<|0j. 
The  boundary  at  |0|  = |0j  resembles  a classical 
spinodal  line.  This  leaves  open  the  possibility  of 
a region  I < 82<  02  of  metastable  gas  or  liquid 
which  behaves  as  an  ideal  single-phase  fluid,  and 
whose  thermodynamic  properties  in  this  region 
continue  to  obey  Eqs.  (B3)-(B6)  rather  tha.i  Eqs. 
(B8),  (B9),  etc. 

In  the  limit  6-6*  (Appendix  BC),  the  boundaries 
0m  and  0„  merge  into  a single  boundary  of  meta- 
stable behavior  along  which  both  the  specific  heat 
and  the  susceptibility  remain  finite,  with  an  in- 
finite derivative  in  the  susceptibility7  and  a cusp 
in  the  chemical  potential  at  this  boundary.  For 
6 <6*  (as  is  the  case  for  the  values  of  6,  ,3,  and  y 
found  in  this  investigation)  the  boundary  |0„|  is 


reached  before  |0J.  For  6->6*,  as  suggested  in 
some  previous  investigations,1*  the  boundary 
|0|  = |0„|  is  not  reached  before  |0j,  and  there  is  no 
accessible  spinodal  line  or  boundary  of  meta- 
stable behavior  in  the  above  sense;  rather,  the 
susceptibility  then  remains  finite  for  all  |0|<|0j, 
and  decreases  to  zero  at  |0|  = |0j,  a situation  that 
seems  unphysical.  It  is  to  be  noted,  finally,  that 
for  the  Ho-Litster9  “cubic -model”  modification 
(end  of  Appendix  C)  that  leads  to  x ^ r~r  as  an  exact 
result,32  the  susceptibility  has  an  infinite  derivative 
and  the  chemical  potential  has  a cusp  just  as  in  the 
restricted  linear  model  in  the  limit  6 — 6 *.25 

APPENDIX  C:  DERIVATION  OF  PARAMETRIC 
REPRESENTATION  FROM  SUSCEPTIBILITY 
DATA 

A.  Natural  variables  fur  use  close  to  critical  point 

The  particular  parametric  representation  postu- 
lated for  use  in  the  analysis  presented  in  this  paper 
can  also  be  derived  from  experiment.  Key  fea- 
tures of  the  experimental  data  are  shown  in  Fig.  2. 
The  coexistence  curve  and  contours  of  constant 
susceptibility  shown  in  Fig.  2 are  not  perfectly 
symmetric  about  the  critical  density.  However, 
to  the  extent  that  the  contours  are  symmetric  (the 
smaller  susceptibility  contours  in  particular  are 
quite  symmetric)  they  can  be  characterized  ap- 
proximately by  the  following  empirical  relation- 
ships11: 

X*(Tp-Tc)-y,  (Cla) 

| Ap|oc  (_aT)b  coexistence  boundary,  (Clb) 

|Ap|oc  (Tt  - T )“(TP  - T,.)6’’  susceptibility  contour. 

(Clc) 

HereAp=p-pc,  AT  -T -Te;  Tp  is  the  maximum 
temperature  on  a susceptibility  contour;  and" 
y = 1.2,  q=  0.5,  ,3  = 0.35. 

For  q =i,  and  with  the  substitution  of  parametric 
variables  r and  9 defined  through  the  relations 
( Tp  - r)(Tp  - Tc)~'  -b282  and  Tp  - Tc-Tir , expres- 


sions  (Cla)  and  (Clc)  become 

X = X(t)  ~Xtr~r 

(C2a) 

and 

A p(r,  0)  =p,rfl  0, 

(C2b) 

with 

AT(r,  0)  = r,r(  1 - 6202). 

(C2c) 

Here  62  is  a dimensionless  constant  and  x, , p,. 
and  Tt  are  constants  with  dimensions  of  suscepti- 
bility. density,  and  temperature,  respectively. 
The  latter  two  equations  [(C2b)  and  (C2c)J  are 
equivalent  to  Eqs.  (1)  of  the  text. 
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It  is  evident  that  in  this  representation,  illus- 
trated in  Fig.  3,  the  parameter  r is  constant  on 
each  contour  of  the  type  (Clc).  The  magnitude  of 
y grows  in  proportion  to  Tp  - Te ; that  is,  r grows 
in  proportion  to  the  “remoteness”  of  the  contour 
from  the  critical  point.  The  parameter  9 is  pro- 
portional to  the  macroscopic  order  parameter  Ap 
at  constant  r;  the  average  “order.”  9.  increases 
from  left  to  right  along  each  contour  of  given  r. 

A constant  value  of  |0|  with  b292>  1 corresponds 
to  the  boundary.  (Clb),  of  the  region  of  coexisting 
phases. 

The  parameter  r as  specifically  defined  above 
is  closely  related  to  both  the  susceptibility  at 
equilibrium  of  the  macroscopic  thermodynamic 
state  and  the  microscopic  range  £ of  transient 
spontaneous  ordering.  According  to  the  approxi- 
mate representations  (Clc)  and  (C2)  of  the  experi- 
mental data  for  C02,  contours  of  constant  r as 
defined  above  are  contours  of  almost  constant  sus- 
ceptibility. In  accordance  with  the  expected  de- 
pendence of  the  susceptibility  on  5,2’3  it  may  be 
conjectured  that  these  contours  are  also  contours 
of  nearly  constant 

It  is  known  that  along  the  critical  isochore 
(Ap  =0)  the  dependence  of  •?  on  temperature  is  of 
the  form3'1'1,26  tj<*( T -Tc)~v.  iThe  exponent  v is 
found  both  theoretically  and  experimentally  to 
satisfy  the  relation  y = 2v  or,  more  exactly,3 
y - (2  - t))v,  with  t) -O.l3,14,26,33).  Consequently,  the 
parameter  r as  defined  above  can  be  given  the 
following  alternative  interpretations: 

Macroscopic:  (x/x,)*l/r,  (C3a) 

Microscopic:  (C3b) 

The  empirical  statements  (Cl)  and  the  assign- 
ment q = i are  approximations.  Consequently, 
not  all  of  Eqs.  (C2)  and  the  deductions  based  on 
them  should  be  expected  to  exactly  characterize 
a real  fluid,  such  as  carbon  dioxide.  Equations 
(Cl)  and  their  consequences,  (C2 )etseq.,  serve, 
rather,  as  simple  approximately  correct  charac- 
terizations of  the  critical  region  that  can  be  per- 
fected through  comparisons  with  more  exact  ex- 
perimental data  and  refinements  of  theory. 

By  using  the  linear  model  to  fit  the  data  we  have 
effectively  retained  equations  (C2b)  and  (C2c)  as 
exactly  correct  equations  of  transformation  from 
conventional  p,  T to  parametric  r,  9 variables. 
However,  we  have  chosen  the  constants  in  those 
transformation  equations  to  meet  our  requirement 
that  the  susceptibility  be  mainly  a function  of  r 
and  thus  maintain  the  approximate  validity  of 
(C2a)  and  of  the  interpretations  based  on  Eq.  (C2a) 
that  have  been  discussed  above. 


B.  Equation  of  state  close  to  critical  point 

Isothermal  integration  of  (C2a),  using  (C2b) 
and  (C2c),  leads  to  a nonanalytic  component  of 
the  chemical  potential  (E2)  (Appendix  E)  which 
includes  as  a factor  an  infinite  series  in  9.u  A 
good  approximate  representation  of  the  result  of 
this  integration,  for  a choice  of  constants  that 
leads  to  6,.,,  *-l.  = +1 . is 

a p(p,T)  = p.(p.  T)-p(pc,r)  = p,»->*8(e-e3), 

(C4) 

with  p,  =const.  (Higher -order  terms  in  the  infinite 
series  in  0 vanish  identically  for  y + ,1  = 1.5.  and 
they  are  still  quite  small  for  1 .5  < y + ,3<  1 .6,  an 
inequality  that  is  almost  certainly  satisfied  for 

co2.) 

The  simple  expression  (C4)  is  equivalent  to  the 
“linear-model”  expression  (Blc).  A second  in- 
tegration of  (C2a)  will  be  found  to  give  an  expres- 
sion for  the  free  energy  (E3)  that,  with  an  ap- 
propriate choice  for  f(pc,T).  reduces  to  a sum  of 
just  three  terms  equivalent  to  Eq.  (B3)  in  the 
approximation  (C4). 

The  simplified  expression  on  the  right-hand  side 
in  (C4)  is  of  the  form  introduced  n priori  by  Scho- 
field in  his  theory  of  the  critical  region  of  magnets 
and  fluids6  (Appendix  B).  Equation  (C2c)  with  a 
general  choice  of  b2  was  also  postulated  by  Scho- 
field. and  Eq.  (C2b)  is  equivalent  to  the  “linear - 
model”  specialization  of  his  general  parametric 
expression  for  the  density  of  a fluid.  The  above 
discussion  shows  that  for  a special  choice  of  62, 
(C4)  |or  equivalently  (Blc)]  can  be  derived  for 
fluids  from  empirical  data,  summarized  by  the 
approximate  expressions  (Cl),  in  a manner  com- 
patible with  retention  of  all  of  expressions  (C2)  and 
(C3)  as  approximately  correct  statements. 

The  use  of  (C4)  as  an  exact  model  equation  re- 
quires a modification  of  at  least  one  of  the  three 
expressions  (C2).  However,  since  the  required 
modification  is  small  for  appropriate  choices  of 
the  constants  of  the  transformation,  we  have  used 
the  simple  "linear  model”  of  Schofield  with  con- 
stants chosen  to  preserve  (C2a)  as  an  approxi- 
mate, but  not  exact,  equality,  in  the  expectation 
that  any  further  modification  of  Eqs.  (C2)-(C4) 
for  use  close  to  the  critical  point  of  C02  will  be 
small. 

An  example  of  a small  modification  to  the  para- 
metric equations  we  have  presented  is  the  “cubic - 
model”  parametrization  of  Ho  and  Litster5  in  which 
a factor  (1  + ctf2)  is  included  on  the  right-hand  side 
of  (C2b)  [or  equivalently  (la)].  With  this  inclu- 
sion, and  with  the  temperature  and  chemical  po- 
tential given  by  (C2c)  and  (C4),  respectively,  Eq. 
(C2a)  can  be  retained  as  an  exact  result  compati- 
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ble  with  the  parametrization  of  Ap,  AT,  and  Ap, 
provided  that32  h2  = b2M  = 3/(3  — 2/i)  and  c = [2(y  +»3)-  3j/ 
(3  - 2ii).  For  the  values  of  y and  J used  in  this 
paper,  6®,  =1.30,  which  is  also  the  value  of  b2  that 
has  been  used  in  this  paper  because  it  provides  the 
best  over -all  fit  to  the  coexistence  curve  and  to 
the  susceptibility  contours  (see  Sec.  II),  and 
c ^0.05. 


APPENDIX  D:  ALTERNATIVE  REPRESENTATIONS  OF 
MEASURED  SPECIFIC  HEAT 

The  result  of  measurements20  of  the  specific 
heat  of  C02  over  the  range  of  temperatures 
4xl0'5<lAT/Te|<10“2  is  summarized  in  Eq.  (8) 
in  the  text.  Apart  from  a factor  T/Tc  =*  1 , that  ex- 
pression is  of  the  form 


£x 


= (A0(-“o+B0)- 


JL 


(Dl) 


where  f = |A7yrc|.  Equation  (Dl)  closely  resembles 
the  theoretical  expressions  Eqs.  (13)  and  (B9)  of 
the  text,  with  the  addition  of  background  contribu- 
tions. A comparison  with  the  experimental  results 
(8)  (using  R =8.314  jr'mor’,  Tc  = 304  K,  Pc 
= 72.9  atm  =7.38  J/cm3,  wpc=0. 468  g/cm3,  and  a 
molar  mass  equal  to  44  g)  leads  to  the  following 
values  for  the  constants  A0,  B0,  and 

A0+  =20.3,  A0_  =38.1, 

Bot  = -12.6,  B0.  = -0.087,  (D2) 

a0  =0.125. 

The  following  modified  form  of  Eq.  (Dl)  has  been 
used  to  represent  the  experimental  data  approxi- 
mately, in  a form  more  readily  compared  with 
theory  when  a = a'  * «0: 


C , /■*  T 

■y  = (A'f-" 


where 


/ - — / " nn 

A’ 


' r 


B'  =A0/2-“o 


■fr“ 


-AT2a 


■B„ 


(D3) 

(D4a) 
(D4b) 

It  is  to  be  noted  that  expression  (D3)  is  equiva- 
lent to  (Dl)  at  all  temperatures  for  a'  = «0  and  that 
it  is  also  equal  numerically  to  (Dl)  for  any  a'  * or0 
at  the  two  temperatures  I = f,  and  t = t2.  Because 
of  these  features,  (D3)  can  provide  nearly  as  good 
a representation  of  the  data  as  expressions  (Dl)or 
(8)  for  appreciable  ranges  of  «'  and  t.  In  particu- 
lar, for  the  choices  f,  =7*  10"5  and  t2  =2. 5x10"’, 
the  agreement  between  (Dl)  and  (D3)  is  ~l<i  for 
values  of  «'  as  small  as  0.08  and  for  temperatures 
in  the  range  4*  10'5<<<  5*  10'\  This  range  of 
temperatures  includes  most  of  the  specific -heat 


data  reliable  to  ~l<k.20  Values  of  A'  and  B'  cal- 
culated from  experiment  for  this  choice  of  values 
for  /,  and  t2  and  representative  choices  of  a'  are 
tabulated  in  Table  V. 

The  corresponding  values  of  A±  and  B t for  use 
in  Eq.  (8)  are  smaller  than  the  A'  and  B'  values 
listed  in  the  table  by  a factor  VPJRTC  = 0.275. 

APPENDIX  E:  SEPARATION  OF  ANOMALIES  FROM 
BACKGROUND 

A rigorous  discussion  of  critical-point  anomalies 
requires  consideration  of  nondivergent  and  analytic 
contributions  to  thermodynamic  properties,  in 
addition  to  the  divergent  and  nonanalytic  ones 
represented  by  the  equations  of  the  parametric 
models.  Detailed  discussions  of  extensions  of 
the  model  equations  have  been  presented  else- 
where.18,19,25  For  present  purposes,  it  is  suffi- 
cient to  imagine  each  thermodynamic  function 
divided  into  two  parts:  (a)  a portion  that  contains 
the  anomalous  rapidly  varying  behavior  typical 
of  the  critical  point;  and  (b)  background  terms 
that  are  small  and/or  slowly  varying  near  the 
critical  point. 

It  is  convenient  to  begin  this  discussion  by  ex- 
pressing the  inverse  of  the  susceptibility  (A7d)  as 


■or1). 


(El) 


Here  (x"’)a  contains  the  essential  singularity 
characteristic  of  the  anomalous  behavior  of  the 
compressibility  in  the  vicinity  of  the  critical  point, 
and  ( x~')b  is  a background  term  that  is  negligibly 
small  at  the  critical  point  and  expresses  all  de- 
partures from  ideal  critical  behavior  that  appear 
at  appreciable  distances  from  the  critical  point. 
Using  this  decomposition,  it  is  possible  to  discuss 
anomalies  that  occur  close  to  the  critical  point 
without  reference  to  thermodynamic  behavior 
farther  from  the  critical  point.  In  this  representa- 
tion the  latter  behavior  appears  explicitly  as  a 
small  perturbation  to  be  added  to  "ideal”  critical - 
point  anomalies;  the  "background”  terms  con- 
tribute only  secondarily  to  the  critical -point  phe- 
nomena treated  in  the  present  investigation,  but 
are  required  to  assure  a smooth  and  orderly  tran- 
sition from  the  universal  laws  of  behavior  charac- 

TABLE  V.  Values  of  A'  and  B'  calculated  from  exper- 
iment for/,  7x  10'5  and  t2  =2.5  v 10'3,  and  for  represen- 
tative choices  of  o' . 


a' 

A\ 

A'. 

B', 

0.08 

45.3 

84.9 

-42.8 

-56.6 

0.10 

31.0 

58.1 

-26.0 

-25.3 

0.125 

20.3 

38.1 

-12.6 

-0.087 
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teristic  of  all  substances  close  to  their  critical 
points  to  the  less  easily  characterized  and  sub- 
stance-specific laws  of  behavior  that  apply  to 
liquids  and  gases  far  from  their  critical  points. 

Successive  integrations  of  at  constant  tem- 
perature lead  to  expressions  for  the  critical  and 
background  components  of  the  chemical  potential 
and  free  energy: 


P = f \~'dp 

= J (X~‘)adp+f  (x~')bdp  + p(pc,T) 

Pc  Pc 

= Pa+Pb’  (E2) 

/=  f p dp  = f padp + f pbdp+/(pc,T ) 

■>  J,’c  ->pc 

= /.  + /».  (E3) 


and 

(E5b) 

APPENDIX  F:  HEIGHT  PROFILES 

In  a fluid  at  thermodynamic  equilibrium  in  the 
presence  of  an  external  field  of  force  (e.g.,  gravi- 
ty) many  o!  the  thermodynamic  properties  vary 
with  location  in  the  fluid.  For  example,  in  the 
absence  of  forces  other  than  gravity  the  pressure 
varies  with  height  z in  accordance  with  the  rela- 
tion dP  = -mpg0dz.  Here  g0  is  the  acceleration  of 
gravity  and  mp  the  mass  density  of  the  fluid.  As 
a consequence,  the  density,  free  energy,  chemical 
potential,  compressibilities,  etc.,  are  also  func- 
tions of  elevation  within  the  fluid.  From  the  Gibbs - 
Duhem  relation,  (A6),  and  the  definition  (El)  it 
follows  that  (a),  lor  constant  temperature  , 35 


The  decomposition  into/,,  /t,  and  p„,  pb  is  made 
in  accordance  with  the  requirements  that  /,  fa, 
and  pa  be  analytic  except  at  the  critical  point 
and  and  fa  correspond  to  (x~l)0. 

Since  p(p,T ) must  be  analytic  for  all  p*pc  (ex- 
cept along  the  boundary  of  the  region  of  coexisting 
phases)  it  follows  that  p(pc,  T)  is  analytic  for  the 
choice  Ap(p,  T)  given  in  (C4).  Consequently,  the 
model  expression  for  Ap.  given  in  (C4)  can  be 
equated  to  p„  in  (E2).  The  term  pt  in  (E2)  then 
includes  the  (analytic)  portion  of  p(p,  T)  close  to 
the  critical  point  that  is  not  dependent  on  density, 
namely,  p(pc,T),  in  addition  to  all  corrections  to 
expression  (C4)  that  are  required  at  temperature 
and  density  distances  remote  from  the  critical 
point. 

Successive  differentiations  of  (E3)  with  respect 
to  temperature  give  expressions  for  the  entropy 
(A 7b)  and  specific  heat  (A7a).  It  is  to  be  noticed 
that  the  nonanalyticity  in  x"1  leads,  via  (E2)  and 
(E3),  to  the  appearance  of  a nonanalytic  component 
of  the  specific  heat;  the  specific  heat  has  an 
anomaly  as  a consequence  of  (and  calculable  from) 
the  anomaly  in  the  equation  of  state. 

Other  thermodynamic  quantities  can  be  divided 
similarly  into  singular  and  background  components. 
As  an  example,  we  cite  the  pressure,  given  by 


— =dp  = -mgBdz  =x~'dp,  (FI) 

whence 

A z-z-zc  = — f dp  - ^ ~ ^ 

>nHo  Jpc 

(F2a) 

where  zc  is  the  “critical  height”  in  the  fluid,  i.e., 
the  height  at  which  p =pc.  If  p(pc,  T)  does  not  con- 
tribute to  pa , then,  from  (E2),  the  integral  on  the 
right-hand  side  of  (F2a)  is  equal  to 

I X"'dp=Ma+  (x~')bdp.  (F2b) 

Jf>c 

If  (b)  the  temperature  is  not  constant,  it  follows 
from  the  Gibbs-Duhem  relation  (A6)  and  the  iden- 
tity 

that 

dP=-mp«°dz--p(^)TdP+(^)dT 

and  hence 


(F3) 

(F4) 


P=pp-/=R0+Ps, 

where 


(E4) 


dz  - 


I / X~'  dp 

+(mpg0)-'(3P/3r),  r. 


(F5) 


(E5a)  where  T,=dT/dz. 


I a ~ PPa  “ J a 
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